Evolution of curvature and anisotropy near a nonsingular bounce 
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We consider bouncing cosmologies in which an ekpyrotic contraction phase with w 3> 1 is followed 
by a bouncing phase with w < — 1 that violates the null energy condition. The bouncing phase, 
induced by ghost condensation, is designed to produce a classically nonsingular bounce at a finite 
value of the scale factor. We show that the initial curvature and anisotropy, though diluted during 
the ekpyrotic phase, grow back exponentially during the bouncing phase. Moreover, curvature 
perturbations and anisotropy are generated by quantum fluctuations during the ekpyrotic phase. In 
the bouncing phase, however, an adiabatic curvature perturbation grows to dominate and contributes 
a blue spectrum that spoils the scale-invariance. Meanwhile, a scalar shear perturbation grows 
nonlinear and creates an overwhelming anisotropy that disrupts the nonsingular bounce altogether. 
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I. INTRODUCTION 

Cosmological models are expected to explain the large 
scale properties of the early Universe: the homogeneity, 
flatness and isotropy of the background, and the nearly 
scale-invariant spectrum of the primordial perturbations. 
These conditions can be achieved in either a rapidly ex- 
panding (inflationary) phase right after the big bang [IH3] 
or a slow contracting (ekpyrotic) phase preceding the big 
bang [11 [5]. In order to have the latter scenario, however, 
the Universe has to undergo a bounce from the contract- 
ing phase to an expanding phase. This transition imposes 
an extra condition on cosmological models, namely, the 
desired properties of the large scale structure created in 
the contracting phase must be propagated safely through 
the bounce. The bounce may be induced by quantum 
gravity effects when the scale factor shrinks to the Planck 
scale near a classical singularity [6 - 8 . Although theories 
of such a quantum bounce (also referred to as a "singular" 
bounce) are not fully developed, it is conjectured that 
the homogeneity and the scale-invariance conditions gen- 
erated before the bounce should pass smoothly through 
the bounce [BH5]. The intuition is that the low energy 
physics that generates structure on large length scales 
should decouple from the physics at high energies and 
short lengths that is responsible for the bounce. 

An alternative approach that has been pursued in re- 
cent years is to realize the bounce in a classically non- 
singular way [151 [TTJ Q25] • In this approach the Universe 
stops contraction and reverses to expansion at a finite 
value of the scale factor a, when classical gravity and ef- 
fective field theories are still valid. The intended advan- 
tage is that it becomes possible to describe with known 
theories the entire evolution of the cosmological back- 
ground and the conditions that emerge from the bounce. 
To have a nonsingular bounce requires violation of the 
null energy condition (NEC). A commonly used approach 
is to introduce a scalar field that undergoes ghost conden- 
sation [15] . The period during which the NEC is violated 
will be referred to as the bouncing phase. 



Some earlier studies [19] of nonsingular bounces have 
not included an ekpyrotic phase. In those cases, the con- 
tracting universe is unstable to the growth of anisotropies 
unless one assumes extreme fine-tuning of initial condi- 
tions; it is also difficult to produce a sufficiently wide 
spectrum of scale-invariant fluctuations. Moreover, some 
earlier treatments [17j [19] of the bouncing phase only 
consider the last stage when the Hubble parameter H 
varies linearly with time and passes through zero. This 
period does not cover the whole bouncing phase during 
which NEC is violated, leaving out important problems 
that arise in earlier stages of the bouncing phase. 

In this paper, we consider nonsingular bouncing cos- 
mologies that have an ekpyrotic phase prior to the bounc- 
ing phase. The ekpyrotic phase stably smooths and flat- 
tens the Universe, and leads to the generation of nearly 
scale-invariant perturbations. After the ekpyrotic phase, 
ghost condensation is triggered to bring the Universe into 
a bouncing phase with w < — 1 that violates the NEC. 
We consider the entire bouncing phase from the start 
of NEC violation at H = to the nonsingular bounce 
at H = 0. We examine perturbations around the back- 
ground evolution and identify problems that break the 
homogeneity, isotropy and scale-invariance conditions. 

In a previous paper [20], we argued that the bounc- 
ing phase strongly modifies the spectrum of curvature 
perturbations. An adiabatic mode, though decaying and 
negligible during the ekpyrotic phase, was shown to grow 
exponentially just before the bouncing phase when the 
equation of state w crosses —1. By the time the Uni- 
verse bounces, the total curvature perturbation becomes 
dominated by this adiabatic mode, which unfortunately 
carries a blue spectrum that is inconsistent with cos- 
mological observations. This result does not depend on 
how the scale-invariant mode is generated in the ekpy- 
rotic phase. To further elaborate this idea, we will show 
that our conclusion does not change with different gauge 
choices, hence posing a serious physical problem for non- 
singular bouncing models. 

In addition, we will demonstrate a new effect that sig- 
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nificantly amplifies the anisotropy. This amplification oc- 
curs gradually through the entire bouncing phase, and is 
not noticed in previous analyses that focused only on the 
last stage of the bounce when H increases linearly with 
time. In our analysis, we find that the scale factor has to 
decrease exponentially during the entire bouncing phase, 
causing the anisotropy to grow substantially before the 
bounce. Consequently, the homogeneity and isotropy of 
the Universe achieved in the preceding ekpyrotic phase 
are lost. Here we consider both anisotropies that have 
classical and quantum origins. The classical anisotropy, 
present at the beginning of the ekpyrotic phase, generally 
comes to dominate over the ghost condensate field before 
the bounce, unless it is tuned to a sufficiently small value 
even before the ekpyrotic phase begins. A more serious 
issue is the anisotropy induced by quantum fluctuations, 
whose amplitude is linked to the amplitude of curvature 
perturbations. This quantum generated anisotropy is 
negligible during the ekpyrotic phase, but becomes over- 
whelmingly large in the bouncing phase. As a result, the 
Universe is inevitably led to chaotic mixmaster behav- 
ior and contracts to a singularity, a — > 0; that is, the 
nonsingular bounce is totally disrupted. 

In Sec.|TlJ we present a sample model that incorporates 
the generic features of an ekpyrotic phase and a nonsin- 
gular bouncing phase. The bouncing phase is described 
in detail in Sec. |III| Then in Sec. |IV[ we demonstrate 
the dramatic growth of curvature perturbation that al- 
ters the power spectrum, presenting our analysis in both 
the comoving and the synchronous gauges. In Sec. [VJ 
we analyze the evolution of the anisotropy, showing how 
it is induced by the curvature perturbation and how it 
disrupts the nonsingular bounce. 
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we argue 



In Sec. 

that both problems of curvature and anisotropy can be 
attributed to the exponential difference between the low 
energy scale of ghost condensation and the high energy 
scale of the ekpyrotic phase. 



II. THE EKPYROTIC NONSINGULAR 
BOUNCING MODEL 

In this Section, we introduce a generic nonsingular 
bouncing cosmology that we will use for our study. Pre- 
vious analyses of nonsingular bouncing models can be 
divided into those with an ekpyrotic contraction phase 
with equation of state w > 1 and those in which the 
contraction phase has w strictly < 1. Without an ekpy- 
rotic phase, the contraction is unstable to the growth of 
anisotropy and the onset of chaotic mixmaster oscilla- 
tions [S]- To counter the instability, extreme fine-tuning 
of initial conditions must be imposed, or the unstable 
contracting phase must be made so short that it is in- 
sufficient to explain the large scale uniformity and den- 
sity perturbations observed today. For these reasons, the 
only cases of interest are those with a long ekpyrotic con- 
traction phase with w > 1, followed by a transition to a 
bouncing phase with w < — 1. 




FIG. 1: The kinetic term P(X) versus X = \{d(f>f in the 
effective Lagrangian for the scalar field cj>. During the ekpy- 
rotic phase, X is in the linear region X ^> X c and increases 
by a factor e 2N . In the transition to the bouncing phase, X 
decreases by an even greater factor to reach X — X c . During 
the bouncing phase, X further decreases to X < X c where 
NEC is violated. 



The transition is hard to achieve by combining two in- 
dependent components with w > 1 and w < — 1 respec- 
tively, because if the w > 1 component dominates during 
the ekpyrotic phase, then the w < —1 component would 
grow much more slowly and never come to dominate dur- 
ing the bouncing phase. A more economical approach is 
to have a single component whose equation of state w 
switches from > I to < — I between the ekpyrotic and 
the bouncing phase. This scenario is achieved naturally 
by having a scalar field that rolls down a steep negative 
potential during the ekpyrotic phase, and then undergoes 
ghost condensation after the ekpyrotic phase is over. In 
this way the entire contracting and bouncing history of 
the Universe is driven by the same field. 

This framework is described, for instance, in the new 
ekpyrotic model [TB], where the effective Lagrangian for 
the scalar field 6 is 



c = y/=g[p(x)-v(<t>)], x = -\{d<f>f 



(1) 



The background metric g^ v is taken to be flat Friedmann- 
Robertson- Walker, 



ds 2 



-di 2 + a(t) 2 (dx 2 + dy 2 + dz 2 ), 



(2) 



and we use reduced Planck units 8ttG = 1. The kinetic 
term P(X), shown in Fig. [lj is canonical, P{X) w X, 
for large A; but has a minimum at a low energy scale 
X c where ghost condensation takes place. The potential 
V(4>) is shown in Fig. [2] where, from right to left, it is first 

approximated by a negative exponential, — Vo e~^ 2 ^ p ^ 
with p< 1; then bottoms out and undergoes a steep rise. 
Unlike in [TBI HI] where V is designed to bend sharply 
at a fine-tuned value above zero, we choose to let it rise 
smoothly as in Fig. [2] The Universe evolves through the 
ekpyrotic phase (m > 1) to the bouncing phase (w < 
— 1), with a transient kinetic energy dominated phase 
(w w 1) in between, as indicated in the figures. 
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FIG. 2: The potential term V(<p) in the effective Lagrangian 
for the scalar field <j). The ekpyrotic phase corresponds to the 
exponential decline from Vek-beg to Vek-end near the minimum 
of the potential. The kinetic phase refers to the quick rise 
from V e k-end to V c ~ 3p | Vek-end | ■ The nonsingular bouncing 
phase occurs at V > V c . 
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where N measures the total number of e-folds of modes 
that exit the horizon during the ekpyrotic phase [13 . 

To have a nonsingular bounce, H has to pass from 
negative to zero, requiring a bouncing phase during which 
H > 0. According to the second Friedmann equation, 



H 



(11) 



then, the nonsingular bouncing phase must have w < — 1, 
a violation of NEC. This is not possible if the kinetic term 
P(X) is linear in X as in the ekpyrotic phase, because 
for the Lagrangian the above equation becomes (Ap- 
pendix pl2j)) 



The ekpyrotic phase starts at an intermediate energy 
scale X ek _b cg when the scalar field rolls down the po- 
tential from Vek-beg) and ends when the field reaches the 
bottom of the potential Vek-end- The negative exponential 
form of the potential leads to a homogeneous attractor 
solution, 

a = a(-t) p , a = a ck _ cnd (-ff Gk _ cnd /p) p , (3) 

(4) 

(5) 



a 



(-*) ' 

^ = v/2p log (v^Sr (-*))' 



where t is negative and increasing towards zero. This 
solution has a constant equation of state, 



w = i - 1 » i . 



(6) 



where typically p ~ 10 2 [2]. There is a scaling relation 
between the kinetic and the potential energies, 



H 2 



P 



1 



3,-1^=^' X =2< 



According to the first Friedmann equation, we have 



H 2 = 



1 



3k 

a 2 + a 6 ' a 3a+™) 



'■'0 



(7) 



(8) 



where k = 0, ±1 represents the spatial curvature, the a 2 
term represents the anisotropy (c/. Section [v]), and the 
dots may include matter, radiation, etc. The ratios of the 
scalar field energy density to other components, including 
curvature and anisotropy, scale as 



H = -XP 



.x 



(12) 



which implies H = -\<i> 2 < for P(X) = X. The idea 
behind the form of P(X) in Fig. [l] is that X exits the 
linear region after the ekpyrotic phase, and enters the 
ghost condensate region X < X c where P.x < 0. To be 
consistent with the ekpyrotic phase, we must have the 
relation X c < X ck _i, cg , which implies 



X, 



X c ^ ^ck-beg ^_ 2N 



ck-cnd 



x, 



ck-cnd 



(13) 



This exponential factor lies at the core of the problems 
that we will show in this paper. 

The transition from the ekpyrotic phase to the bounc- 
ing phase is mediated by a brief kinetic phase between 
the two. As in Fig. [2] after the ekpyrotic phase, the po- 
tential V(4>) rises sharply in order to slow down the field 
(f> and reduce X from X e k- G nd to X c . Accordingly, the 
equation of state changes from w ^> 1 to w = — 1. This 
transient phase lasts much shorter than a Hubble time, 
during which the total energy is almost conserved. Let 
V c be the value of the potential when X reaches X c , then 
we have 



Vc = 3H 2 



3H, 



ck-cnd 



3p\V e 



ck-cnd | 



(14) 



which means V c is smaller but of a similar order of mag- 
nitude to |Kk-end|- 

The bouncing phase begins when X further decreases 
to less than X c , whereby the equation of state w falls 
below —1. The kinetic energy of the scalar field, 



T{X) = 2XP X - P 



(15) 



-3(l+u>) 



Pk.cr, 



cx (-fT 2 CX H 2 



(9) 



Therefore, after the ekpyrotic phase, all other compo- 
nents including the initial curvature and anisotropy are 



becomes negative, as shown in Fig. [3] The bounce occurs 
when X decreases to the point where the negative kinetic 
energy T cancels the positive potential energy V , giving 
H 2 = |(T + V) = 0. Here the negative kinetic energy 
does not incur ghost instability around X w X c as long 
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FIG. 3: The kinetic energy T(X) = 2XP X - P for the scalar 
field (j). It is a linear function of X in the region X ^> X c , but 
becomes negative at X < X c . 



remain negligible. However, the bouncing phase cannot 
be made arbitrarily short. By the end of the ekpyrotic 
phase, the Hubble parameter H is negative and exponen- 
tially large, so there must be a long period with H > 
in order for H to increase to zero. 

To determine how long the bouncing phase really lasts, 
we solve the equation of motion for the scalar held, 



T x X + 6HP tX X + V 4> (/) = 0. 



(18) 



Using (16 1, a general solution to Eq. (18) can be found 



that describes a nonsingular bounce. First notice that 
Eq. (16) implies, upon integration from X c , 



|XP X |>P, forA<A c . 
Therefore the kinetic energy becomes 

T{X) = 2XP^ X - P « 2XP X ■ 



(19) 



(20) 



The potential energy and its gradient are nearly constant 
during the bouncing phase, 



v{4>) « v c , v^v 4 



(21) 



as T x — P,x +2XP xx i s positive [TS]. This is satisfied 
if P{X) grows quickly as X decreases; in particular, if 



XP 



,xx 



> -P. 



X 



for X < X c 



(16) 



which implies a large slope in T. In that case the bounce 
would happen at an X very close to X c where the ghost 
condensation model is valid. 

We note that the speed of sound, 



P^x 
Tx 



P 



x 



2XP 



XX 



p 



X 



(17) 



takes the canonical value 1 in the ekpyrotic phase when 
P{X) is linear, but becomes negative in the bouncing 
phase when X < X c . Condition (16) is equivalent to 
\Cg\ <C 1 for X < X c , which limits the rate at which 
gravitational instability develops when NEC is violated. 
In our analysis, we will avoid the gravitational instabil- 
ity by assuming an extremely small value of \c 2 s \ in the 
bouncing phase, as quantified in Section |TV} Then we can 
focus on problems of growing curvature and anisotropy 
that arise regardless of the gravitational instability. 



III. THE BOUNCING SOLUTION 

Before studying the evolution of perturbations, a more 
detailed description of the bouncing phase is needed. One 
major issue is the curvature and anisotropy terms in the 
Friedmann equation (JsJ) . The curvature and anisotropy 
are exponentially suppressed during the ekpyrotic phase, 
but they start to grow faster than the scalar field energy 
when w < — 1 during the bouncing phase. It is commonly 
assumed that the bouncing phase can be made as short as 
a few Hubble times, so that the curvature and anisotropy 



< Nt 



-N 



and similarly for V a,, where Aibp is the duration of the 



bouncing phase as found below. The factors 



and 



| y 1 * 1 | are taken to be much greater than 1 yet much 
less than e*; otherwise one has to fine-tune the steep- 
ness of the potential to super Planckian scales, which is 
unphysical as we shall argue in the end of the paper. 
We are now in a position to solve the equation 



T x X + 3HT + V^4> = 0. 



(23) 



The solution can be described in three stages accord- 
ing to whether the friction or the gradient term dom- 
inates. At the very beginning of the bouncing phase, 
\H\ ps \H C \ = ^/Vc/3 and \T\ is small, hence the fric- 
tion term is negligible. Next, the negative kinetic energy 
\T\ increases and the friction term overtakes the gradi- 
ent. Finally, very close to the bounce we have T ~ — V c 
but | if | becomes small, so that the friction term is again 
subdominant to the gradient term. Among these three 
stages, the first and the last are both very short and do 
not contain interesting features. Indeed, since the fric- 
tion and the gradient terms have the same sign during 
the bouncing phase, we have 



Denote 




^) y/2X c V c /3 



(24) 

(25) 
(26) 
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then the first stage goes from T = until T ~ Tb p _beg, 
which lasts less than a Hubble time, 



th 



p-bog 



tr 



{-V,<t> c W 2X c/3V c 




hence the scale factor does not change much, ab p _beg 
«c ~ a c k-ond- Similarly, the third stage begins from T 

^bp-end 

the bounce can be bounded by 



3-ffbp-end , and the time it takes to reach 



p-cnd 



3#bp-cnd < {-V^/V C ) 2 QX C 
f-V. A \ I 2X 



v c Jy P x ek . end \h c \ ~ \H, 



-N 



(28) 



which is much less than a Hubble time, implying w 
abp-ond- Finite factors like p and ( ~y* c ) are neglected 
in these estimates. 

Therefore the bouncing phase mainly consists of the 
middle stage during which the friction term dominates 
and the gradient is negligible. Under this condition, the 
equation of motion simplifies to 



f + 3HT = 0, with H 2 = i(T+ V c ). 
The bouncing solution is readily given by 



cosh 2 (||i7 c |(i- t )) ' 
H=\H C \ tanh(§|# c |(i-i )) , 



(29) 



(30) 
(31) 



where the bounce is formally at t = to. This solution 
applies between the moments ibp-beg and ibp-end, which 
correspond to \T\ ~ |T bp -bog| and \H\ ~ |i?b P -end| respec- 
tively. From the expression (25) and use Eq. (30), we 
have 



COSh 2 (§|if c |(*bp-beg-*o)) ~ 

which gives 

^0 ^bp-beg 



-N 



X, 



ck-cnd 



N 

3Th7 



(32) 



(33) 



Similarly, from the expression (26) and use Eq. (31), we 
find 



tanh (§|-ff c |(*bp-eiid-io)) ~ 
which gives 

^0 ^bp-end 



X, 



-N 



X, 



ck-cnd 



2e 



3H C 



(34) 



(35) 



Therefore the middle stage, hence the whole bouncing 
phase, lasts approximately for a period 



A^bp ~ ^bp-end ^bp-beg 



N 1 



(36) 



This duration eliminates any hope to complete the 
bounce within just a few Hubble times. As a result, 
we expect anisotropies to grow significantly during the 
bouncing phase. From the bouncing solution (31), the 
scale factor a(t) scales as 

acxcosh 2/3 (§|i? c |(t-i )) cx |T|~ 1/3 . (37) 

Before the end of the bouncing phase, it contracts by a 
factor 



^bp-end 
^bp-beg 



-^bp-end 
Tbp-beg 



\ v c J \x ck _ c 



< 



e 3 



N 



(38) 

Therefore, naively, the anisotropy term in the Friedmann 
equation ([8| would increase by a factor 



( (j2 )bp-end 
(cr 2 ) bp -bcg 



(^bp-end \ 
O-hp-b™ ' 



x r 



> 



2N 



(39) 



which completely cancels the suppression ( |10| it has 
experienced during the ekpyrotic phase. That is, the 
anisotropy has returned! 

To be more precise, during the ekpyrotic phase the 
suppression of anisotropy happens in such a way that 
the anisotropy term itself remains almost constant, while 
the scalar field energy increases by an exponential factor 
e 2N . In the bouncing phase, however, the anisotropy 
term grows by the same factor e 2N , while the scalar field 
energy decreases rapidly towards 0. Thus there must 
be a point near the end of the bouncing phase where 
the anisotropy contribution to the Friedmann equation 
^ overtakes the scalar field energy. After this point, 
the nonsingular bouncing solution (31 ) that assumes the 
dominance of the scalar field becomes invalid. To see 
how the exponential growth of anisotropy may disrupt 
the nonsingular bounce, we need to study perturbations 
around the background evolution found above. 



IV. CURVATURE PERTURBATION AND 
POWER SPECTRUM 

Consider linear perturbations about the metric |2]), 

ds 2 = a(r) 2 [ - (1 + 2A)dr 2 + 2(B ti + S l )drdx i (40) 
+ ((1 - 2ip)5ij + 2E 4j + 2F {h3) + 2h ij )dx i dx j , 

where r is the conformal time, dt = a dr. We will denote 
conformal time derivative as ' = 4z = aj^, and spatial 
Laplacian as V 2 = d l di. The functions A, B, ip and E 
represent the scalar perturbations; Si and Fj, with S\ — 
F\ = 0, represent the vector perturbations; and hy, with 
h\ = h 1 -^ = 0, represent the tensor perturbations. The 
definition of these perturbative quantities are subject to 
gauge transformations of the coordinates, as discussed in 
Appendix IB] The constant time hypersurface is given by 
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the normal — (— a(l + A), 0). The perturbation of its 



intrinsic curvature is (Appendix (B5)) 



(41) 



hence is referred to as the curvature perturbation. The 



shear of the hypersurface is given by (Appendix (B7|) 

^=«[(^"I^W)+ ( 7^. ) ], (42) 
where the scalar and vector contributions are 



E' -B , 
F'i - Si , 



(43) 
(44) 



which will be related to the anisotropy. The evolution of 
these perturbative quantities is determined by the Ein- 
stein equation and the equation of motion for the scalar 
field 4> and its perturbation 5<j>, as shown in Appendix [B| 
The evolution of the curvature perturbation involves 
only the scalar perturbations. To agree with obser- 
vations, the contracting phase of the Universe has to 
generate a nearly scale-invariant spectrum of curvature 
perturbations. Several mechanisms have been proposed 
that can give rise to such scale-invariant curvature per- 
turbation, including the multi-field cntropic mechanism 
[H)l |2"T] and the single-field adiabatic ekpyrotic mecha- 
nism [22l [23] • Once the generated scale-invariant mode 
exits the horizon, it is expected to be conserved on su- 
perhorizon scales. Meanwhile, there exists an adiabatic 
contribution to the curvature perturbation that is time- 
varying according to the equation of state. This adiabatic 
contribution generally has a blue spectrum, and is sub- 
dominant to the constant scale-invariant modes on large 
scales. However, a serious problem arises when the equa- 
tion of state w abruptly changes in the transition from 
the ekpyrotic phase to the bouncing phase. It happens 
that an adiabatic mode grows exponentially as w drops 
past —1, and surpasses the scale- invariant mode even on 
large scales. 

To observe this problem, consider the curvature per- 
turbation in the comoving gauge defined by Sip = 0, as 
is done in |20j . The curvature perturbation in this gauge 
is described by the gauge-invariant quantity (Appendix 
dB2§) 



K = ip + H^ . 



(45) 



The Fourier mode lZ k with comoving wavenumber k 
obeys the equation (Appendix (B33)) 



n'l + 2-n' k + c 2 s k 2 n k = o, 



(46) 



where z = ay —2H /c 2 H 2 . After the mode exits the hori- 
zon, the k 2 term can be treated perturbatively, and the 
equation is formally solved by expanding in orders of k 2 , 



n k = nf-k 2 



dr 



dr c 2 z 2 lZ k 



(47) 



The leading order 1Z k is the solution to the equation 
without the k 2 term, which contains two general solu- 
tions, 



K k a) = C 1 (k)+C 2 (k) 



dr 

^2 



-^const 



Kf. (48) 



These two terms are the leading adiabatic contributions 
to the curvature perturbation; the fc-dependence of the 
dominant term determines the spectral index of adiabatic 
perturbations on large scales. 

The fc-dependence of the constants C\ , C 2 can be found 
by matching to initial conditions in Minkowski vacuum 
when the mode originated from quantum fluctuations 
deep inside the horizon. Introducing the canonical vari- 
able v k = z lZ k , Eq. ( 46 1 becomes 



H = o. 



(49) 



z 

During the ekpyrotic phase, c 2 — 1 and the scaling so- 
lution (|7|) gives z (— r) Tzr ? , so the freeze-out horizon 



scale is given by 



P(l-2j>) 

(1 -p) 2 T 2 



(50) 



Then the solution to Eq. (49) is given by Hankel func- 
tions, 

v k (r) = ^[AHi 1 \x) + BH^{x% (51) 

where x = k(—r), and v = | — jzr^- At early times 
when the mode is deep inside the horizon, x 3> 1, the 
solution should be matched to the vacuum solution for a 
Minkowski background [24 , 



1 



-ikr 



2k 



(52) 



Using the asymptotic behavior of Hankel functions, 
~2 



,±i(x- 



f) 



x — > 00, 



(53) 



we can fix the constants B = 0, A — \fjj: + ^ ■ 
Therefore, neglecting the constant phase factor, we have 



(54) 



After the mode exits the horizon, the solution approaches 
the other limit, a;< 1, where it allows the expansion 



Vk 



.2 u T{v) 



( 



.cos(vk)T(— v) 
2 v ir 



(55) 
x -> 0. 
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From this, the expression for lZk becomes 

/ / vl-J- 

iirp \—T) 2 1 -p 



+ i 



8 (1 — p) 1 -p a l ~*> 
cos(z^7r)r(l — v) 



2 V TTV 



— ^f-k 2 + + v \ k 2 ^(-t) 1 -?, 
2a 1 -p 2a 1 -p 



(56) 



in the limit of small p, but keeping the exact powers. 
These leading terms correspond to the solution ( 48 ) , 



+ TZ^ = C 1 (k) + C 2 (k) 



o ai{-2H) 

-2 

pa 1 -? 



(57) 



Ci(fc) + C 2 (fc)^— ((l-p)(-r)) 1 - 



2-6p 

Comparing the coefficients, one obtains 

2a 1 -f 
i 



C 2 (fc) 



k 2 i- 



(58) 
(59) 



Since fc 3 / 2 |d(fc)| - fc, and fc 3 / 2 |C 2 (fc)| - fc 2 , both 72.£ onst 
and 72,™* terms have blue spectral indices. 

To produce a scale-invariant contribution to the cur- 
vature perturbation, we may, for instance, invoke the 
entropic mechanism described in [16] . In this mecha- 
nism, an independent scale-invariant contribution 72,| c " lnv 
is generated from entropic perturbations by an additional 
scalar field at the end of the ekpyrotic phase, where 



72.sc-mv 



f 



k 3/2 



(60) 



Thus the total curvature perturbation is the sum of the 
terms, 



Kl ot 



^rjconst 



(61) 



By the end of the ekpyrotic phase, according to 
Eq. |57l 



I ck-cnd 



k I ek-cnd 



cm, 

C 2 {k)p 2 
2(a 3 H) 



ek- 



C 2 (k) 

^ a ek-end 



X, 



ck-cnd 



(62) 
(63) 



Hence 72' fe nt is suppressed relative to 72.™ nst by a factor 



7ef 




p k 


•^const 


ck-cnd 


(aH) c k-cnd 



x k 



-N k 



X, 



ck-cnd 



where Xk is the kinetic energy at horizon crossing and 
N^, is the remaining number of e-folds of the ekpyrotic 
phase after the /c-mode exits the horizon. Due to its blue 



spectrum, the 72.™ nst term is in turn dominated by the 
scale-invariant term, 



V {~r)- V 


■y^const 




k 






•^sc-inv 


ek-cnd 


\Q>H) ck-cnd 





Xi. 



Hence the integral term is sub-subdominant, 



-N k 



(65) 







k 


2 


x k 


•^sc-inv 


ck-cnd 


(ai?)ck-cnd 




-^ck-cnd 



(66) 

Therefore, the total curvature perturbation right after 
the ekpyrotic phase is 



72* fc ot 



ck-cnd 



^sc-mv 



f 



fc 3/2 



(67) 



which confirms that the entropic mechanism produces a 
scale-invariant spectrum. We note that the lZ s k c ~ lnv and 
72.™ nst terms remain constant afterwards, but the inte- 
gral term 72.' fe nt is time-varying even on large scales. The 
total curvature perturbation TZ k is conserved and scale- 
invariant on superhorizon scales only if the time- varying 
piece remains negligible. 

We now demonstrate that, in fact, the 72.jf* term grows 
rapidly in the kinetic phase when H increases from neg- 
ative to and w crosses —1. Indeed, from Eqs. (12 1 and 



(17), the integral term can be written as 
2 H 2 



TZi nt 



C 2 {k) 



-2H 



dt = C 2 {k) 



1 H 2 dX 



2XT X X 



(68) 

As the field passes the bottom of the potential and climbs 
up the other side, its kinetic energy quickly decreases to 
X s» X c where P(X) becomes nonlinear. During this 
rapid phase the friction term in Eq. ( 18 ) can be neglected, 
hence 



72-jf 



C 2 (fc) 

C 2 (k) 
C 2 (k) 



1 H 2 



dX 



a- 



2Xy 

2 



3« 3 



1 /Hi 

2a 3 V^ c 



/2X 
x " dX 
V2X3 

1 

V2X~c ' 



(69) 



where the integral in the first line is dominated by con- 
tributions from near the upper limit, and so we approxi- 
mated the slowly varying quantities by their values there. 
Compared to the value in ( 63 ) at the end of the ekpy- 



rotic phase, the integral term has now been exponentially 
amplified, 



l^-fc | ck-cnd 



x, 



k-cnd 



-WV X 

It exceeds the scale-invariant term by a ratio 
Xk I X, 

r^ 7 



(70) 



^sc-inv 



k-cnd 



ck-cnd > e N-2N k > ^ 

x c ~ 

(71) 
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for a wide range of modes with Nk < N/2. 

Thus, modes that exit the horizon in the second half of 
the ekpyrotic phase, including all modes within current 
observable horizon, are overwhelmed by the integral term 
which carries a blue spectrum, 



(72) 



in contradiction to observations. Moreover, assuming 
that the amplitude of the scale-invariant modes in ( 71 ) 



matches the amplitude measured by the cosmic mi 

j£ 

2.2 x 1CT 9 

tude A 2 , » e 2N ~ m *A 2 n 



crowave background (CMB), A 2 ^ 



|^sc-i„v|2 



, then the blue modes will have an ampli- 
nv , which is already nonlinear 
for modes with Nk < N/2 — 5. Note that these modes re- 
main outside the horizon after the ekpyrotic phase, since 
the horizon continues to shrink in the kinetic phase. This 
can be seen by writing 



where we used the fact that \H\ < \H C \ in the bounc- 
ing phase, and assumed that AX <C X c in order for 
the ghost condensate model to be valid. Hence the total 
change of the integral term is much less than the value 
( 69 1 obtained in the kinetic phase right before the bounc- 



ing phase begins. 

Our analysis differs in an important way from [171 119) , 
which considered similar models and concluded that the 
comoving curvature perturbation changes negligibly near 
the bounce. In those cases, the bounce is analyzed in the 
limit that the Hubble parameter varies approximately 
linearly with time. From our bouncing solution (31 1, this 
corresponds to the period when \t — to\ <C l/\H c \^so that 



H « %H 2 (t — t c 



H 



H 



(78) 



This is within the last e-fold of the boun cing phase, where 



<C 1 and a w ab p _end given in Eq. ( 38 ) . During this 



linear regime the curvature perturbation changes by 



a 2 H l 



,(^)_ eysfT^) 

-2+y(l + W )-^(l + W ) 2 



which is dominated by the first term as w 



-1. 



(%) 



( a ^)ek-end- 



(73) 



(74) 



Since we assumed a steep potential in the kinetic phase, 
this expression is generally larger than the value at the 



end of the ekpyrotic phase, 
Eq. |50|. 



\z"/z\ a (aH)i k _ end /p, from 



It remains to show that the dominantly blue curvature 
perturbation is maintained through the bouncing phase 
without significant changes. By Eq. (68), the change in 
the integral term is 



AK k nt 



C 2 (k) 



H 2 dX 
2a 3 T ~X~ 



(75) 



We observe that during the whole bouncing phase we 
have the relation 



H 



< const 



1 



3o 3 

ou ck-cnd 



V, 



-v. 



1 



(76) 



AH k nt 



C 2 (k) 

c 2 (fc) 

a 3 

"bp-end 

C 2 (fc) 



bp-end 



-2H 



^dt 



(79) 



J 3 -H 2 (t-t ) 2 dt 



a ek-end 

C 2 (k) 
a 3 

"ck-end 



-v. 



x, 



ck-end 



x r 



\-V. A JJ2} 



Thus, as found by [T71 [TU] . 1Z changes very little from 
its value (691 at the beginning of the bouncing phase. 



Indeed, we can simply see from the expression ( 79 1 that 
the integral term is decaying as \H\ — > at this last stage 
of the bouncing phase. However, we emphasize that the 
comoving curvature perturbation grows large at a much 
earlier stage just before the bouncing phase begins. This 
period was ignored in those previous studies, and so the 
problem with the blue spectrum was missed. 

Finally, let us take care of the gravitational instability 
in the bouncing phase when c 2 < 0. During the bouncing 
phase, perturbation modes may reenter the horizon. For 
example, in the case of a slowly varying c 2 (i.e. c s /Hc s <C 
1) consistent with Eq. (16), the freeze-out horizon scale 
corresponding to the bouncing solution (31) and (37) is 
given by 



The equality holds for the bouncing solution due to 
Eq. ( 29 1 and the scaling relation ( 37 1 ; the value of the 



constant is found at the beginning of the bouncing phase 
from Eq. (25). The inequality is true in the short stages 
before ibp-beg and after £bp 



,_ond as a result of Eq. (24) 
Therefore we can estimate Eq. ( 75 ) as 



ATli 




(77) 



a 2 H 2 



1 1 

I ~ 2cosh 2 (§|ff c |(i-i )) 
9 

2sinh 2 (|| J ff c |(i-i )) 



(80) 



which goes from negative to positive and crosses zero at 
to — t s» 1A/\H C \. Therefore all modes briefly reenter the 
horizon around this time, and exit again when t — to — > 
approaching the bounce. According to Eq. ( 49 1 , modes 



inside the horizon grow unstable instead of undergoing 
oscillations. Since the growth rate is proportional to 



9 



|c s |fc, this instability could be tamed if \c%\ is small and 
the duration inside the horizon is short. The danger- 
ous modes are those on small scales, which reenter the 
horizon early in the bouncing phase and stay inside until 
very close to the bounce. Their duration inside the hori- 
zon can approach the upper bound ( 36 1 , or in conformal 
time, 



IfcAr 



(81) 



which may create a problem for the modes with Nk < 
N/3. To avoid this gravitational instability, the speed of 
sound has to be 



< 



N 



(82) 



This is the same condition under which the leading term 
( 48 1 dominates the expansion in Eq. ( 47 1 and our com- 



putations follow. We note that even when gravitational 
instability is suppressed under the condition ( 82 ) , the 



problem with the exponential growth of the comoving 
curvature perturbation still occurs just before c 2 s becomes 
negative. 

Having traced the evolution of the comoving curvature 
perturbation through each stage of the contracting phase, 
let us now comment on the curvature perturbation in 
other gauges. Under a general coordinate transformation 
x 1 - 1 — > x^ + £ ^ , the curvature perturbation ip transforms 



as (Appendix (B22|) 



ip -> tp + aHC 



(83) 



An immediate consequence is that the curvature pertur- 
bation is gauge-invariant at a nonsingular bounce where 
H = 0; hence results computed in different gauges must 
agree as they approach the bounce. Therefore the expo- 
nential growth of curvature perturbation that we found 
in the comoving gauge is truly physical and should exist 
in other gauges as well. However, although the curvature 
perturbation in different gauges end up with the same 
value at the bounce, their patterns of growth can be quite 
different during the contracting phase before the bounce. 
For example, in the longitudinal gauge (Appendix [B]) 
the Newtonian potential $ grows exponentially during 
the ekpyrotic phase, but undergoes no abrupt change 
in the transitional kinetic phase, then grows further in 
the bouncing phase to reach the same value as TZ at the 
bounce. In the above computations, we chose to study 
the comoving curvature perturbation TZ because it will 
be conserved outside the horizon in the expanding phase 
including the reheating period, and therefore it is directly 
related to the density and temperature fluctuations that 
are observed today. It is also a convenient variable to 
study because it is almost conserved during both the 
ekpyrotic and the bouncing phases, whereas the expo- 
nential growth can be clearly confined and ascribed to 
the transition between the two phases, during which the 
equation of state w changes rapidly. 

Since the comoving curvature perturbation TZ becomes 
exponentially large right before the bouncing phase, an- 



other natural concern is whether the perturbative compu- 
tation breaks down at this point. We address this ques- 
tion by noting that the perturbation theory is valid if 
there exists one particular gauge in which the pertur- 
bations of all relevant physical quantities remain small; 
then the variables in other gauges can be formally de- 
fined and related to the quantities in this particular gauge 
through gauge transformations. In our case, the syn- 
chronous gauge plays this role. It can be readily checked 
that the physical quantities, including the curvature per- 
turbation -0 S , the shear perturbation a s (see Section [V)), 
and the matter perturbations Sp s , SP S , (p + P)Su s , all 
remain finite during the transitional kinetic phase in this 
gauge. In particular, the curvature perturbation ip s in 
the synchronous gauge can be related to TZ through (Ap- 
pendix dB49l)) 



ip s = TZ-H 



TZ 
H 



dt' 



(84) 



After the k mode ip sk exits horizon, we may insert the 
leading terms of TZ/, from ( 48 1 to find (for brevity we will 
omit the subscript k from here on) 

ip s « d(k) + C 2 {k) [ dt'H f dt"^L . (85) 
Jo Jq 2a- i H 

Due to the presence of H in the outer integral, the sec- 
ond term remains small when H — > 0, so that the cur- 
vature perturbation is well behaved in this gauge when 
the equation of state w crosses —1. Hence it is indeed 
legitimate to carry our computations beyond this point 
into the bouncing phase. 

Nevertheless, even the synchronous curvature pertur- 
bation ip s has to grow exponentially near the bounce to 
match the value of TZ. To prove this, let us evaluate 



Eq. (85) in the bouncing phase as follows. First denote 



+ ^ nt = Ci(fc) + C 2 (k) / dt'H ■ I(t') , (86) 



where I it') is the inner integral in (85). In the kinetic 
phase, similar to Eq. (|69|, 1(f) can be computed as 



(87) 



Here this huge growth as t' — > t c is tempered by H — > 
in the outer integral in ( 86 1 , hence does not show in ip 1 ^ 




Then in the bouncing phase, as for Eqs. (75) and (77), 
I(if) can be found to stay nearly constant, 



AI 



x 



H dX 



< 





f v \ 




AX 


7T~3~ 

■^ck-end 




\H c \y/2X c 


2X C 



«I(*c). 
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Therefore it simply contributes a constant factor to the 



outer integral in (86), which leads to 



4's 



Ca(fc) / dt'H 



C 2 {k) 



3 a 



ek-i 



\-V.d,JJ2. 




Clearly when H — > near the bounce, ^' nt approaches 
the same value as R lnt given by Eq. (69), which is ex- 



ponentially large and nonlinear for a wide range of k 
modes. From the above derivation we also see that ij) s 
becomes exponentially large only when \H\ <C \H C \ near 
the bounce. This corresponds to \t — to| <C l/|-ff c | m the 
bouncing solution (31), which is within the last e-fold 
before the bounce. We shall see below that during this 
time the anisotropy grows large and nonlinear, too. 



V. ANISOTROPY AND NONSINGULAR 
BOUNCE 

Anisotropy is described by the shear of the constant 
time hypersurface. In Appendix [X] we derive the gen- 
eralized local Friedmann equations, and show that the 
anisotropy term exactly corresponds to the squared mag- 
nitude of the shear. Here we illustrate a simpler case with 
a flat, homogeneous, but anisotropic metric, 

ds 2 = -N 2 dr 2 + a 2 %(dx i + (3 l dT){dxi + (3 j dr), (90) 

where the lapse A/", the shift a 2 (i l , and the spatial metric 
a 2 "fij depend only on r. We choose the scale factor a(r) 
such that det7ij = 1. The extrinsic curvature of the 
constant time hypersurface with normal n M = (—A/", ) 
is given by 

Kij = rij-i = |^7ij + a.^ , (91) 

where it is decomposed into the expansion 9 and the shear 
Oij by using the projection tensor 7^ = g^ v + n^n v , 

fl _ i 



''>:.) 



;k 



(92) 
(93) 



and % 



is the conformal Hubble parameter. The 



squared magnitude of the shear is given by 



a 2 = \a»ai 



(94) 



For a perfect fluid at rest, e.g. a homogeneous scalar 
field, the Einstein equations can be written as 



n — n° — 1 n' 2 nr 2 
P - -(* - - ° J 

3p = G', = ~ff -9 2 -3a 2 



= G\ 



1 si f^ih 
3°J U k 



1 j< 



6a\ 



(95) 
(96) 
(97) 



In physical time dt = Afdr, the Hubble parameter is 
H = H/Af = 9/3, and the first two equations reduce to 
the Friedmann equations, 



H 2 



H 



-\{p + p)-a 2 



The third Einstein equation becomes 



a 



3Ha\ 



0, 



j 1 3 
which implies that the shear scales as 

1 2 1 



(98) 
(99) 

(100) 
(101) 



Thus the anisotropy term in the Friedmann equation Q 
is precisely the squared shear of the constant time hyper- 
surface. More general cases with nonflat and inhomoge- 
neous metrics are dicussed in Appendix [A] 

For our analysis of the bouncing dynamics, we 
again consider the perturbatively inhomogeneous and 
anisotropic metric d40|, 



ds 2 = 



- (1 + 2A)dr 2 + 2(B i + S i )drdx i 



((1 - W)5i 



2E, 



2F, 



(i,3) 



2hij)dx % dx 3 



Small anisotropies can be analyzed by studying the shear 
perturbation of the constant time hypersurfaces. The 
normal to these hypersurfaces is n M = (— a(l + A),Q), 
and the shear is given by (Appendix (|B7[)) 



B 



(E'-B)) 



(id)) 



(102) 



It involves both the scalar and the vector perturbations, 
which evolve independently at linear order. The scalar 
part of the shear is coupled to the curvature perturbation 
through Einstein equations, while the vector part scales 
in a simple manner with the scale factor a. Therefore, we 
expect the anisotropy to grow significantly whenever the 
curvature perturbation becomes large, or when the scale 
factor shrinks exponentially as in the bouncing phase de- 
scribed in Section |III| In that case, the rise of a large 
anisotropy would substantially alter the dynamics of the 
contracting Universe and potentially demolish the non- 
singular bounce. 

Let us first look at the vector perturbations that fol- 
lows very simple behavior. The vector part of the shear 
perturbation is given by 



- a ( F (iJ) ~ S (hi)) = acr (M) 



where 



r v _ 



(103) 



(104) 



and (T Vj i = 0. The variable is gauge-invariant by 
itself, and evolves according to the equation (Appendix 



2H a) 



0. 



(105) 
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which implies 



and hence 



Vi 1 v 1 v 1 



(106) 



(107) 



in agreement with (101). Unlike in an expanding uni- 
verse, the vector perturbation of , if nonzero, grows in 
the contracting universe as the scale factor a decreases. 
Therefore, in general, the vector perturbation should not 
be overlooked in analyzing bouncing cosmologies [25], 
since it naturally creates anisotropy and may disrupt the 
bounce. As we argued in Section |III| the exponential 
growth of anisotropy in the bouncing phase may cause 
the breakdown of perturbation theory near the bounce. 
However, since the vector perturbation is not sourced by 
the scalar field, it will remain negligible if it is extremely 
small initially. In this sense, the problem with the vec- 
tor part of the anisotropy can be resolved by fine-tuning 
initial conditions. 

Suppose that the universe is homogeneous but there 
exists some initial anisotropy <7^ k _ bcg at the beginning of 
the ekpyrotic phase. We would like to determine if, af- 
ter exponential amplification in the bouncing phase, it 
may dominate the energy density and prevent the non- 
singular bounce. For this purpose, we shall include the 
anisotropy term in the Friedmann equations and reana- 
lyze the bouncing behavior of the Universe. Under the 



approximations (20) and (21), the Friedmann equations 
become 



3H A = P<p + cj z «T + V c +<t a , 
-T 



H 



-XP, X - a 2 



(108) 
(109) 



In order to reach a nonsingular bounce in the presence 
of a nonzero a 2 , p^ m ust decrease to a negative value to 
cancel the a 1 in (108). Meanwhile, H must stay positive 
before H hits zero; hence, the kinetic energy T has to 
be sufficientl y ne gative so that its magnitude is larger 
than 2cr 2 in (109). From Eq. (37) we know that in the 
bouncing phase T scales as 1 /aF, slower than a 2 which 
grows as 1/a 6 . Therefore if a 2 starts with too large an 
initial value, it may overtake the scalar field energy and 
prevent the bounce. 

Here we derive a general condition for the nonsingular 
bounce to happen. Our only assumption is that X, and 
hence T, is monotonically decreasing with time in the 
bouncing phase, which is guaranteed by the steep poten- 
tial that slows down the field (f>. Let the anisotropy be 
<j\ at some time t\ when it is still much smaller than the 
scalar field energy p\. Introduce the variable 



(110) 



which starts from q = 1 and increases as a decreases. 
Since a 2 oc a~ 6 oc q 2 , we have 



dq 



2q. 



For the scalar field, from Eq. ( 23 1 , we have 
fa + 3HT = 0, 

and hence 

dp4, _ T Ti 
dq q ~ q ' 



(111) 



(112) 



(113) 



since T is monotone decreasing. Therefore, Eq. ( 108 ) 
satisfies 



d(3H 2 ) < _\T 1 \ 



2a 2 q, 



dq q 
which integrates to give 

3H 2 < 3H 2 - \Ti\ \ogq + a\{q 2 - 1) 



(114) 



(115) 



To have a bounce at finite a, or finite q, a sufficient condi- 
tion is that the RHS has a root in q > 1, or equivalcntly, 
the minimum of the RHS is less than 0. The minimum 
is at q 2 = \Ti\/2a 2 , where the RHS equals 



Pi 



2 



2 g 2a 2 



Requiring this to be less than 0, we find the condition 
I Ti| __, 



o\<^ exp{- 



\Ti\ 



!}• 



(116) 



This condition applies for any point t\ in the bouncing 
phase, in particular the last stage where Eq. (29) and 
hence T oc 1/a 3 does not hold. Nonetheless, since it is 
sufficient at any point t\ , we may find the mildest condi- 
tion by a proper choice of t\ . Let t\ be within the middle 
stage where we can use T oc 1/a 3 . Extrapolating back to 
the beginning of the bouncing phase, we need to have 



^bp-beg 



< 



a 



exp{- 



bp-beg 



204-1711) 
\Tx\ 



1} 



IT, 



bp " bcgl rxnh - ^ 
2 i Ti | expjl m 



(117) 



The mildest condition is established at t\ — > tbp-ond where 
\Ti \ ^ |r bp _ond| ~ V c , which imposes 



^"bp-beg 



< 



IT, 



bp-beg 



2U C 



-V 



<t>c 



Vc 



2 Xc 

3e 



(118) 



Therefore, since a 2 is nearly constant through the ekpy- 
rotic phase, the condition on the initial value of Cgk-beg 
is simply 



^■2 < v 

°ck-bcg A c , 



(119) 
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which has to be much less than p<p ck _ hcg ~ 3pX e k-bc g - In 
other words, to prevent the anisotropy from becoming a 
problem during the bouncing phase, the initial conditions 
must be fine-tuned so that the ekpyrotic phase is highly 
isotropic to begin with. 



Condition (119) applies for not only the vector per- 
turbation but any initial anisotropy that scales as 1/a 6 , 
including any initial scalar contributions. However, for 
scalar perturbations, large anisotropy can arise even if 
there is no initial anisotropics classically, because the 
scalar perturbation is continuously sourced by the quan- 
tum fluctuations of the scalar field. We demonstrate 
that the scalar shear perturbation is induced during the 
ekpyrotic phase and grows exponentially large near the 
bounce. 

The scalar part of the shear perturbation is given by 



= a{{E% - B tij ) 



„ , - i*yV 2 (tf -B)) 

2 - s ^ (120) 



where 



t s = E' - B 



Under the coordinate transformation 



transforms as (Appendix ( |B23 )) 



s ^ 



(121) 



(122) 



and thus is not gauge-invariant. In the comoving gauge, 
cr s is described by the gauge-invariant quantity 



(48) plus the expansion (47) of TZ, we find, up to order 



c 2 (k) cm r* 

cm ( a 



(127) 



dr. 



t r 2 H 2 



2a 3 B 



dt' 



aB 



o a 2 Jo 



c 2 B 2 
dt" s 



2a 3 B 



Along the same lines as for Eqs. ( 69 ) and ( 85 ), we see that 
the third term increases exponentially as H — > 0, while 
the fourth and all higher order terms are well behaved 
and finite. Compared to the leading term, a c becomes 
dominated by the third term which surpasses the leading 
term by a factor 



c c ok- 



2fc 2 
3(a#)L, 



( v c x m, 



ek-end ^ N-2N k 

2Xr 



(128) 



Note that this is the same exponential factor as in ( [71] ), 
by which the comoving curvature perturbation TZ grows 
when w — > — 1 for the range of modes with < N/2. 

However, this huge growth of the comoving shear per- 
turbation does not imply the breakdown of perturbation 
theory at this point. As we have seen in Section |IV[ in 
the synchronous gauge the curvature perturbation ip s re- 
mains small throughout the contracting phase until very 
near the bounce. Accordingly, the shear perturbation 
would also remain small in this gauge. The synchronous 
shear perturbation a s can be related to the comoving 
shear perturbation a c through the gauge transformation 
(Appendix |B50|) 



(123) 



1 



n 

B^ 



dt' 



(129) 



The evolution of a c is coupled to the comoving curva- 
ture perturbation TZ through the equations (Appendix 



(B3CKB32)) 



2H<7 C 



TZ!_ 
H 



TZ = 0, 



TZ' + 



w -w 



V 2 {TZ + na c ) = 0. 



(124) 
(125) 



Thereby the scalar shear and curvature perturbations 
source each other to grow from quantum fluctuations to 
classical perturbations when the modes exit the horizon. 
From the second equation above we can solve for a c in 
terms of TZ, 



^TZ'- U 

c 2 k 2 n 2 n 



(126) 



Since we know from Section [IV] that the comoving curva- 
ture perturbation TZ undergoes exponential amplification 
right before the bouncing phase when w crosses —1, we 
expect the comoving shear perturbation <r c to follow the 
same growth. Indeed, using the leading order solution 



After the mode exits the horizon, we use Eqs. (127) and 
((57} to find, up to order Oik 1 / 2 ), 



c 2 (k) cm f 



C 



dt 



c 2 B 2 



a 2 \B J 2a 3 H 



1 ,aB 



n JO 



c 2 B 2 
dt" s 



2a?B 
t 



C 2 (k) 



c 2 B 



dt' 
2a 3 B 



C 2 (k) d(k) r , C 2 (k) 



dt 



,aH 
IP 



11 Jo 

C 2 (k) Ci(fc) 
2 



c 2 B 2 
2a 3 B 



dt' 







H 
IP 



c 2 B 2 \ 
2a 3 B J 



2k 2 a 2 
C 2 (k) 



dt'a 

a z Jo 

t ft TT ft 

dt'aB / dt"-^ 
Jo H Jo 



c 2 B 2 
dt'"^^. (130) 
2a 3 B 



From the first to the second line it is clear that the grow- 
ing term in a c is absorbed by the gauge transformation 
term, leaving behind a n in tegral that is well behaved as 
H ->• 0, similar to Eq. H. 
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Moreover, using similar calculations that lead to (89 1, 
<7 S in the bouncing phase can be evaluated as 



C 2 (k) Ci(fc) 
2k 2 a 2 a 2 
C 2 (k) 



dt't 



Jo HMZ\-V 4> )JT) 



C 2 (k) 
2k 2 a 2 



Ci(fc) 



-"■ek-end 

C»(k) (-v c /v t4>c 



3n 3 

° a ck-cnd 



dt'a 



C 2 (k) C\(k)a hp _ hcg 1 C 2 {k) (-K/V^J 
a 2 \H C \ a 2 3a 2 k _ cnd \H c \y/2X c ' 

(131) 



2k 2 a 2 



The integral in the second line is found by using 
(371, which rapidly converges to the constant value « 
a bp-bcg/|-ffc|- From the last line we see that all three 
terms scale equally as 1/a 2 near the bounce, which co- 
incides with the vector perturbation ( 106 1 and likewise 
leads to 



(132) 



that resembles Eq. (101) in the homogeneous case. We 



also find that the third term in ( 131 1 dominates over the 
first (and second) term by a ratio 



2k 2 (-VJV^) 
3«c 2 k-e„d \H c \V2X~c 



(aH) 2 



ck-cnd 



x t 



ck- 



,N-2N h 



x r 



(133) 



for modes with Nk < N/2, similar to (128) but only at a 
much later time in the bouncing phase. 

We note that the scalar shear perturbation is generated 
by quantum fluctuations and cannot be fine-tuned away. 
To see whether the resulting anisotropy grows nonlinear, 
we estimate the size of the anisotropy by computing the 
self-correlation function of the scalar shear perturbation, 



{S)a%{x)) 



d 3 k 1 



(2tt) 3 2a 2 
1 



( 



±A- k 



2 ^s 



S |2 ,„6 



67r 2 a 2 



k°dk 



(134) 



The integration is carried over the modes that exit the 
horizon in the second half of the ekpyrotic phase, i.e. 
those that ha ve < Nk < N/2 and are dominated by the 
third term in ( 131 ). Modes with Nk > N/2 are negligible 
since the power spectrum of er s is deeply blue, whereas 
the modes with Nk < remain inside the horizon un- 
til very close to the bounce, contributing only vacuum 
fluctuations that can be eliminated by renormalization. 
Therefore, as we approach the putative bounce, a lower 



estimate of the anisotropy is 



(fff)ok-ond 



<(^) 2 } > 



67r 2 a 6 



C 2 (k) {-VJV^) 



IOStt 2 ^ ai k _ cnd H 2 X c 



3a 2 k _ cnd \H c \V2Xo 

(a-ff)ek-Bnd 



k b dk 



C 2 (k)\ 2 k b dk 



(-V c /V^) 2 a 



IOStt 2 ^ ai k _ cnd H 2 k _ cnd X c 



ck-cnd /„ rr\8 
~~?T~ \ aI1 Jck-cnd 



"ek-end / V c \ 

a 6 \-V.sJ 



Vc. \ 2 ^ck-end 7^4 

j£ ck-cnd ' 



(135) 



where we used the expression (59) for C 2 {k). At t 
ibp-end, we have 



) ) bp-end 



> 



bp-beg 2N tt4 
e ti, 



a 



ck-cnd 



2 AN 



(136) 



bp-end 



In order for the anisotropy not to become a problem, 



condition (116) requires that 



Vc. 



^bp-end < Y e 

This can be satisfied only if V 2 e 4N < V c , or 



Vc. < 



-4AT 



M, 



PI 



(137) 



(138) 



in reduced Planck units. However, in order for the 
scale-invariant modes to match the observed amplitude 
10~ 5 , the potential energy at the end of the ekpy- 



p 



rotic phase needs to satisfy \/2V c \ i _ cnc \/p • 
which implies 



10- 3 M 2 ! 



V c ps 3pV c 



ok -end 



p 2 lO~ b M ] 



pi ■ 



(139) 



For a typical value of p ~ 10~ 2 H], the condition (138) 
is incredibly far from being satisfied. 

Since the scalar shear perturbation er s is not gauge- 
invariant, we must explain our gauge choice for comput- 
ing (135). Recall that condition (137) is obtained for the 



homogeneous case by using the Friedmann equations of 
the Hubble parameter H . In the inhomogeneous case, the 
Friedmann equations should be generalized to equations 
of the local expansion 9 of the constant time hypersurface 
(see (A38 A39) in Appendix [A"|. In general, the expan- 
sion 6 is not uniform on the hypersurface, and its gradient 
along the hypersurface is determined by a supplementary 
equation (A40). In order to apply condition (137), it is 



preferable to adopt the constant mean curvature gauge 
described in Appendix |A| In this gauge Eq. (A40) van- 
ishes identically, so that is indeed uniform and equal to 
3H on the constant time hypersurface. This gauge co- 
incides with the synchronous gauge in the homogeneous 
case. For linear perturbations, this gauge reduces to the 
uniform Hubble gauge defined in Appendix[B] The scalar 
shear perturbation in the uniform Hubble gauge can be 
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computed through (B57l and is found to be equal to the bouncing phase according to Eq. (39 1. This cre- 



the synchronous shear perturbation a s to le ading order, 
which is what we used to compute Eq. (jl35h above. 



It is worth contrasting the uniform Hubble gauge with 
the longitudinal gauge in which the scalar shear pertur- 
bation cr s is set to zero identically. The vanishing shear 
perturbation in the longitudinal gauge does not mean 
that the bouncing phase is free from problems. Accord- 
ing to Eqs. (B51) and (B17l, the variation of the local 



Hubble parameter over the constant time hypersurface is 
given by 



Sn = -{U'-H 2 )8u 



(140) 



In the longitudinal gauge, thoug h a° = 0, the velocity 
perturbation Su — er c (Appendix (B40l) grows exponen- 
tially due to Eq. (128 1. Therefore the local Hubble pa- 



rameter is extremely inhomogeneous over the hypersur- 
face, which means the nonsingular bounce is ill defined 
in this gauge even when the background Hubble param- 
eter H vanishes. In the uniform Hubble gauge, however, 
8% = on every constant time hypersurface, hence the 
bounce would be reached simultaneously at the hypersur- 
face on which H = 9/3 = 0. Therefore the nonsingular 
bounce is well defined in this gauge, and the huge growth 
of anisotropy is indeed physical. 

In a consistent perturbative analysis, the a 2 term in 
the Friedmann equations should be negligible to linear 
order. The exponentially large ((c s ) 2 ) that we found in 
( 135[ ) implies that the perturbation theory is no longer 
valid. Therefore we expect the anisotropy generated by 
quantum fluctuations to become nonlinear and dominate 
over the scalar field energy before the bounce. Once the 
anisotropy takes over, the Universe will be driven to a 
BKL like contraction phase that ends in an extremely 
inhomogeneous and anisotropic singular crunch [551 HZ] ■ 



VI. CONCLUSION 

We have examined a cosmological model of a contract- 
ing universe that smoothly connects an ekpyrotic phase 
with a nonsingular bounce. We have shown how the non- 
singular bouncing phase spoils the flatness, homogene- 
ity, and isotropy condition, as well as the scale-invariant 
perturbations generated during the ekpyrotic phase. We 
have identified four different effects that create problems 
for nonsingular bouncing models: 

• gravitational instability: This problem is caused by 
c 2 becoming negative in the bouncing phase. One solu- 
tion is to have c 2 remain exponentially small during the 
bouncing phase, as in (82 1. Alternatively, it has been 



suggested [TH1 HE] to restrict the bouncing phase to one 
Hubble time; this requires tuning (— V jt j, c /V c ) to be expo- 
nentially large, which is not viable, as shown below. 
• regrowth of initial anisotropy: The initial 
anisotropy suppressed during the ekpyrotic phase is ex- 
ponentially amplified to even greater magnitudes during 



ates a problem for the nonsingular bounce unless the 
initial anisotropy present before the ekpyrotic phase is 
fine-tuned to an extremely small value given by (1191, 
e.g. by having a dark energy dominated expanding phase 
preceding the ekpyrotic phase, as in [TOl 111] . 



• blue spectrum of curvature perturbations: As 

first discovered in |20j , a sub-subdominant component of 
adiabatic curvature perturbation generated from quan- 
tum fluctuations in the ekpyrotic phase grows to domi- 
nate over the scale-invariant component as the bouncing 
phase begins, resulting in a blue power spectrum that 
is inconsistent with observations. This problem occurs 
when the equation of state w passes through —1, and it 
persists even if the gravitational instability is suppressed. 

• growth of quantum induced anisotropy: The 
scalar shear perturbation generated from quantum fluc- 
tuations in the ekpyrotic phase grows exponentially dur- 
ing the bouncing phase and contributes a large anisotropy 
that dominates the energy density and prevents the non- 
singular bounce. This problem is due to both the equa- 
tion of state w passing through — 1 and the scale factor a 
decreasing exponentially in the bouncing phase. It per- 
sists even if the first two problems are avoided. 

According to Eqs. (70 1, p8]), and (1351, it would ap- 



pear that the last two problems could be resolved if the 
factor (— V^ c /V c ) is made to be at least of order e N ; the 
first two problems would not occur either if this condition 
were achieved, which amounts to keeping the duration of 
the bouncing phase within one Hubble time. How ever , 



this condition requires |V^J > V c e , and by Eq. (139) 



the value of \V^ c | would exponentially exceed the Planck 
scale. It is then unavoidable to consider quantum gravity 
effects in this approach, which defeats the purpose of the 
nonsingular bounce to avoid those effects. 

It is worth noting that none of the above four prob- 
lems appear in the singular bounce, as in the cyclic model 
[TTJl [TT] . The first two problems disappear because in the 
singular bounce c 2 is always positive and the horizon is 
forever shrinking. Moreover, in the singular case, w re- 
mains > 1 throughout the contraction phase and the en- 
ergy density is increasing all the way to the bounce, so the 
curvature and anisotropy never grow to dominate. But 
in the nonsingular approach considered here, w must fall 
below —1 and X must fall below X c , and it is during this 
phase that the curvature and shear perturbations grow 
to make the universe inhomogeneous and anisotropic. 

In the introduction to this paper, we noted a de- 
coupling argument suggesting that the bouncing phase 
should not affect the large scale density perturbations 
or the isotropy of the Universe, as long as the bounce 
occurs at much higher energy densities than during the 
ekpyrotic phase when the perturbations are first gener- 
ated. Our conclusion does not contradict this intuition, 
but instead supports it: the failure of the nonsingular 
bounce considered here is precisely due to the fact that 
the bounce is mediated by ghost condensation, a low en- 
ergy effect rather than a high energy one. 
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Indeed, from Eqs. (70 1 and (39) it is clear that both or in matrix form, 



growths of curvature and anisotropy are proportional to 
the same factor X e k_ eIu j/ X c , the ratio between the high 
energy scale near the end of the ekpyrotic phase and the 
low energy scale associated with the ghost condensation. 
In order to consistently merge the ekpyrotic phase with 
the bouncing phase through ghost condensation, these 
energy scales have to satisfy the hierarchy 



X c « X, 



ck-bc 



< x, 



ck- 



(141) 



Therefore the exponentially large factor is inherent in 
such ekpyrotic nonsingular bouncing models. Our calcu- 
lation illustrates how this same factor naturally shows up 
and leads to the four problems listed above. 

Therefore, we conjecture that any form of nonsingular 
bouncing models would suffer from these problems unless 
the hierarchy ( 141 ) can be relaxed. This may happen 



either by reducing the energy ratio between the beginning 
and the end of the ekpyrotic phase, or by raising the 
energy scale of the bounce. The first possibility is to seek 
a nonlinear realization of the ekpyrotic phase in which X 
varies within a limited range. The other possibility is 
to induce a bounce with physics of a much higher energy 
scale. We are currently attempting to construct examples 
of both types. 
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Appendix A: Generalized Priedmann equations 

1. spatial hypersurfaces and timelike congruences 

Consider a spacetime (A4,g) foliated by a family of 
spacelike hypersurfaces {E Ti t€R}, we may construct 
a coordinate system that is adapted to the foliation as 
follows [55] • Take r to be the natural time coordinate, 
then the lapse function Af is defined as 



Af = ( - (07 



,-1/2 



(Al) 



Let n be the timelike unit vector normal to the constant 
time hypersurfaces, then by definition = (— Af, 0). 
The spatial coordinates are fixed by choosing a shift vec- 
tor /? defined as 



P = d T - Afr, 



(A2) 



which implies 

system, the metric can be written as 
ds 2 = -Af 2 dr 2 + llJ (dx i + P l dr){dx 3 



Under this coordinate 



■ftdr), (A3) 



9»v = 

where j3i 
given by 

in/ 



.9oo 9oj 

9i0 9ij 



-Af 2 +/3 k /3 k ft 
Pi lij 



(A4) 



7ij/3 J . Accordingly, the inverse metric is 



5 oo g0j 

giO gij 



1 

JV -5 

p* 

AT 5 " 



AT 5 



P'P 3 



(A5) 



where ^ is the inverse of the 3-matrix 7^. 

The intrinsic curvature of the hypersurface S is given 
by the Ricci tensor ^Rij associated with the induced 
spatial metric 7^ on the hypersurface. 7^ can be pushed 
forward to form a projection tensor 7 in Ai, 



(A6) 



The extrinsic curvature of the hypersurface £ is then 
given by 

K»v = 7%7 A ^« ; a , (A7) 
where semicolon ; denotes the covariant derivative com- 



patible with the metric ( A3 ) . It can be shown that is 
symmetric, K^ v = K v(11 and tangent to the hypersurface, 
K^ lJ n l/ = 0. The mean curvature of the hypersurface is 
1/3 of the trace K = = K\. 

The 4-dimensional Ricmann tensor R^ Xv k. °f ^ ne s P ace_ 
time A4 can be related to the intrinsic and extrinsic cur- 
vatures of the hypersurface S through the Gauss-Codazzi 
relations [2S], especially, 



K 



= ^R af} + KK aP 

{3) R + K 2 



2R flu n f *n u + R 
R^n^Ya = -K\ a 



K l3 K, 



(A8) 

(A9) 
(A10) 



where | denotes covariant derivative associated with the 
induced spatial metric 7. 

Consider the congruence of the integral curves of the 
timelike normal vector n. Such an integral curve can be 
regarded as the worldline of an Eulerian observer [28] . 
who is defined to have a 4-velocity equal to n. The co- 
variant derivative of the timelike vector n M can be kine- 
matically decomposed as [29] 



where 



a„ = n„ = n n 



fi ; v 



(AH) 



(A12) 



is the acceleration of the Eulerian observer, and the ex- 
pansion tensor 9^ and the vorticity tensor u„ v are sym- 
metric and antisymmetric respectively. It can be shown 
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(e.g. [30] ) that for the unit vector n normal to the hy- 
persurface S, the vorticity tensor uj^ v vanishes, and the 



expansion tensor 



a„iij, is equal to the ex- 



trinsic curvature K^ v . Decomposing 9^ further into the 
trace and the traceless parts, we have 



3" IV 



cr^v - a a n 



where the volume expansion 

= T^< 



(A13) 



(A14) 



is equal to the trace of extrinsic curvature K, and the 
shear is traceless and symmetric. Both the shear 
and the acceleration a M are tangent to the hypersurface, 
a^n v = a u n v = 0. For the general metric (A3) and the 
normal vector = (— A/", 0), the expansion, shear, and 
acceleration can be explicitly computed to be 



1 

w 



7« 7y , k fi k , (A15) 



(A17) 



Using the kinematic decomposition, the Gauss- 
Codazzi relations (A8 A9 A10) can be written as 



= ^R 



a/3 



2R llv n ih n v + R = (3) i? + |6» 2 - 2a 2 



R n^-i" 



a\n 



(A18) 
(A19) 
(A20) 



where a 2 = \<j % i<Jij is the squared magnitude of the shear 

[23 EH ESI- 

To study the kinematic evolution of the timelike con- 
gruence, let us first introduce the Fermi derivative ^ L 
[2"0] that propagates a timelike vector V along its inte- 
gral curve, ^jjV = 0. The Fermi derivative of a vector 
field X with respect to V is defined as 



ds 



X^ = A M — V^V V X V + V»V v X l 



while for a covariant vector u it is 
D F 



ds 



= w„ - VuV v u v + VuV u u) v , 



(A21) 



(A22) 



and similarly for tensors. It reduces to the covariant 
derivative e F e V v X ix . v if the integral curve 

of V is a geodesic, i.e. when = 0. 

The Fermi derivative of the expansion tensor 9^ v with 
respect to n is given by [29] 



Df 9 



Rn\vKn X n K 



9 liX 0\ - 2n^6 v)x a > 



+ 7 M 7 „a(\; K ) + 



(A23) 



The trace of this equation gives the Raychaudhuri equa- 
tion, 



= -R XK n x n K - \9 2 - 2a 2 + a%. (A24) 
And the traceless part of the equation gives 

+ <Vd„ + (R XK n x n K + 2a 2 - a x . x ) . 

(A25) 



Multiplying by 7*"7 I/ - and using Eqs. ( A18 A19 ), we find 



F i 



B^tr*. = 'f'fjRuu - ^R% - ea^ - unyjr 



\Qn l a,j + a l .j + a l aj + d l rij 

i 7 4 ,(i? + R^n"n v - ^R + a" ;/t ). (A26) 



2. local Friedmann equations 

The dynamics of the timelike congruence is determined 
by the Einstein equation relating the spacctime curvature 
and the matter stress-energy tensor. Consider the matter 
source to be a perfect fluid with stress-energy tensor 



T^u = (p+ P)u^u v + Pg^ 



(A27) 



where p and P are the rest energy density and pressure, 
and u is the 4-velocity of the fluid. The fluid velocity 
relative to the Eulerian observer is 



u 

U=- n 



(A28) 



where the Lorentz factor T is given by 

T = -n^u^ = (l-U^y 1/2 . (A29) 

Accordingly [25], the energy density as measured by the 
Eulerian observer is 

E = T ilv <n! i n v = T 2 (p + P)-P , (A30) 

the momentum density as measured by the Eulerian ob- 
server is 



p a = -T^n^ v a = (E + P)U a 



(A31) 



and the stress tensor with respect to the Eulerian ob- 
server is 

S»v = l\l\T XK = (E + P)U^U„ + P 7 ^ . (A32) 

The trace of the stress-energy tensor is given by T = 
= — p + 3P, while the trace of the relative stress 
tensor is 5 ee = T + E. The Einstein equation can 
be written as 



(A33) 
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the trace of which gives 

R = -T = p — 3P 



(A34) 



Projecting the Einstein equation along the normal vector 
n, we find 

R^n" = E + \ {—p + 3P) , (A35) 
whereas the mixed projection with n and 7 gives 

R^j v a =- Pa = -(E + P)U a , (A36) 
and the full projection onto the hypersurface £ gives 



yyyi/ H , = l s < i-§7i(-/'+ 3p ) 

= (£ + J P)ir i [/ J + i 7 > 



(A37) 



Finally, we can use the above equations to write 
Eqs. (|A19|, (tA24b, ||A20| and dA26b as 



{ l9r = l(E-^R + a>), 



3 U \n 



(A38) 
(A39) 
(A40) 



and 



D F . 



9o l 



{E + P)U l Uj - {3) R l 



n l <jj k a k - \0n l aj + a 1 ,j + a l a,j + a 1 
#AE-p-WR + a*). 



3 U J 



(A41) 



Eqs. | |A38| and jA40| ) are equivalent to the ( g) an d (?) 
components of the Einstein equation; Eq. ( A39[ ) cor- 



responds to the trace of the (*■) components, whereas 



Eq. (A41| corresponds to the traceless part. 



These equations closely resemble the Friedmann equa- 
tions for a homogeneous universe. The expansion 9 can 
be considered as 3 times the local Hubble parameter [33] , 
then Eqs. (A38l and (A39) are the local Friedmann equa- 



tions. Note that according to Eq. ( A40 ) the local Hubble 



parameter may vary from point to point on the constant 
time hypersurface. In the homogeneous case, the fluid 
stress-energy tensor 7),„ is diagonal; the expansion 9 and 
the shear cr*. do not depend on spatial coordinates, while 



the acceleration a M vanishes. Accordingly, Eqs. (A38 
A39 1 reduce to the Friedmann equations, 



H 2 



H 



Up-^R + 3 ), (A42) 
-l(p + p) + mR-a 2 . (A43) 



plifies to 



Eq. (|A40|) vanishes identically, whereas Eq. (A41) sim- 
: -SHcj 1 , - PQJZ*. + §<**,• < 3 >H . (A44) 



3. gauge choices 

We may choose certain gauge for the coordinates where 
equations ((AT5WAI7} and (|A38HA4l]) simplify. First let 



lij = a 2 jij 



(A45) 



where the scale factor a is such that det^ = 1. The 
conformal Hubble parameter is therefore T~L = ^. Note 
that here the scale factor a may depend on the spatial 
coordinates, and so does T-L. Then the expansion, shear, 
and acceleration become 



M Ma 



— — i k{l ff] + 



2Afd 



,k > 



Afa 2 



(A46) 

.TV" 

(A47) 
(A48) 



To specify a gauge, the time slicing can be fixed by either 
choosing a foliation condition for the whole spacetime, or 
by starting from an initial Cauchy surface and choosing 
a lapse function Af. Then the spatial coordinates can be 
fixed by choosing a shift vector f3 l . 

The comoving gauge is defined such that the folia- 
tion {E T } is everywhere orthogonal to the fluid velocity 
u, i.e. the Eulerian observer is comoving with the fluid, 
so that n = u. This gauge exists when the fluid does 
not have vorticity, which is true for a scalar field. In this 
gauge the relative velocity U vanishes, so the Eulerian 
observer simply measures the comoving energy density 
and stress tensor, 



E = P , Snu = Pluu 



(A49) 



We may further fix the spatial gauge by choosing j3 % = 0, 
so that the expansion and the shear simplify to 



1 

277 



- ik 



7 ikj 



(A50) 



Here the expansion 9 may depend on spatial coordinates, 
hence the local Hubble parameter is not uniform on the 
constant time hypersurface. Also, since the fluid is in 
general not free streaming, the acceleration a M does not 
vanish. The local Friedmann equations become 



¥ 



2 a i\j 



P) 



3 U ;m 



(A51) 
(A52) 
(A53) 



and the shear a l j evolves according to 



D F , 



9a% - (3) i? 4 , + U\-( (3) i? - a%) (A54) 



u o-jka — ^au dj + a . ■ + a dj + a Uj 
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The synchronous gauge is defined such that Af = 
1 and f3 l = 0, which fixes the coordinates given the 
choice of an initial Cauchy surface and spatial coordi- 
nates therein. In this gauge the acceleration vanishes 
according to Eq. (A17|, so the integral curve of n, or the 



worldline of the Eulerian observer, is a geodesic. Then 
Eq. ( A39 1 becomes 



(A55) 



and Eq. (A41) simplifies to 



a) = -9a) + {E + P) U l U - ( ®& J -Itf^E-p- {3) R) , 

(A56) 

where the Fermi derivative reduces to the covariant 
derivative. 

Another useful gauge is the constant mean curva- 
ture gauge in which the foliation is such that the mean 
curvature 9/3 is constant on each time slice. In other 
words, the Hubble parameter H is uniform on each con- 
stant time hypersurface, and does not depend on the spa- 
tial coordinates. In this gauge, the Friedmann equations 
dA38l |A39b become 



H 2 = l(E-h i3) R + a 2 



and Eq. (A40) becomes a constraint on the shear, 



a\ {j = -(E + P)U t 



Appendix B: Linear perturbations 



(A57) 
(A58) 

(A59) 



On the perturbative level, the general metric ( A3 ) can 



be expanded about the flat FRW metric as, to linear 
order, 



ds 2 



a(r) 2 - (1 + 2A)dr 2 + 2{B ti + S l )drdx l (Bl) 



((1 - 2^)5i 



2E, 



2F (m) + 2h ij )dx i dx j 



where A, B, ip and E represent the scalar perturbations; 
Si and Fi, with S i i — F i i = 0, represent the vector 
perturbations; and hij, with h\ = h % j i = 0, represent 
the tensor perturbations [32]. Comparing to (A3), we 
identify the lapse, shift and the spatial 3-metric to be, 
also to linear order, 

Af = a{l + A), (B2) 

pi = a 2 (B, i + S i ), (B3) 

7y = a 2 ((1 - 2il>)5 ii + 2E tij + 2F {hj) + 2%) . (B4) 

Accordingly, the constant time hypersurface has a time- 
like normal vector 

curvature of the hypersurface is given by 

4 



(-0(1 + A),0). The intrinsic 



&R 



(B5) 



whereas the trace of the extrinsic curvature K, or equally 
the expansion 9, is given by 



1 



9 = - [37Z(1 - A) - 3*1/ + V 2 (E' - B)] . 



(B6) 



Note that the expansion only involves the scalar pertur- 
bations. The shear of the constant time hypersurface is 
given by [33] 

Vii=*{{Kii- B M)-¥^ 2 {E'-B))+a(F' {i>j) -S m ) 



aia , 



i<%vV) 



a a 



v 



where the scalar shear perturbation is defined as 

a s = E' - B , 
and the vector shear perturbation is 

^ = F[ - S t . 

For the matter source, consider a scalar field cj. 
Lagrangian C = P{4>,X), where X — ~^(d(j)) 2 . 
corresponding stress-energy tensor is 



(B7) 



(B8) 



(B9) 

with 
The 



(BIO) 



which takes the form of a perfect fluid (A27l with pres- 
sure P, energy density p = 2XP,x — P , and velocity 
u u = d u 4>/\/2X . In accordance with the homogeneous 



FRW background metric, the scalar field is assumed to 
be homogeneous on the background level as well. There- 
fore the background values 4>, X, and hence p, P, are all 
functions of time r only, and = (— a, 0). Since the 
shear vanishes at this level, the background equations of 
motion are fully captured by the Friedmann equations, 

H 2 = \p=\{2XP, x -P) , (Bll) 

H = -±(p + P) = -XP x . (B12) 

Then let us consider a small perturbation 8((>(t, x ) 
about the background value 4>(t), which also creates a 
perturbation about X(t), SX = 2X(-A + s -&-). Accord- 
ingly, the perturbations in velocity, energy density, and 
pressure are given by, to linear order, 



Su 
5p = 

5P = 



= —aA , Sui = adu .i , Su = ^ — , (B13) 
\ {%' -V. 2 ) \ \ (8v! + USu -A)- 3HSu] , 



(B14) 



~(U'-U 2 ) 



Su' + 2H 



2y 



(H'-n 
{w-n 2 ) 



Su- A 



(B15) 



where the speed of sound is c 2 — , which only needs 
be kept to zeroth order. Here we have assumed <// < to 
agree with our model in the text. 
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The equations of motion for the linear perturbations 
are given by the perturbed Einstein equation, SG^ — 
ST^ . (Alternatively, one may impose the conservation of 
energy and momentum |35j . plus the equation of motion 
for the scalar field.) At linear order, the scalar, vector, 
and tensor perturbations evolve independently and can 
be treated separately. In the following we discuss the 
scalar and vector perturbations. 



1. scalar perturbations 

The equations for the scalar perturbations are 

2 r 



a 

2 r 

O 2 

2 
2 



m(ip'+UA) - V 2 (V> + H(E'-B)) 
(i/Z+nA)] = { P + P)6u , 



tp-A+ (E'-B)' + 2H(E'-B) 



= 0, 



(ip'+HA)' + 2H(i>'+HA) + {U'-H 2 )A 



8p, (B16) 

(B17) 

(B18) 

= SP, 
(B19) 



where the matter perturbations are given in terms of 
the s calar field perturbation 8<f> through Eqs. (B13 B14 
B15[ ). These equations correspond to the scalar modes 
of the (° ), traceless (^), and the trace (\) compo- 
nents of the Einstein equation. There are 4 degrees of 
freedom in these equations, namely A, ip, a s —E' — B, 
and d(f>. However, as we shall see below, these variables 
are not independent since there exists one gauge redun- 
dancy. Therefore, after gauge fixing, we are left with 3 
physical degrees of freedom. Accordingly, given the mat- 
ter perturbations in the form (B13l ~ (B15), one of the 



equations in (B16l ~ (B19) is redundant 



Consider the infinitesimal coordinate transformation 



(B20) 



where for the scalar mode, £j can be written as £j = 
From the transformation rules of the metric tensor 
and the scalar cj>, one finds that the perturbative quanti- 
ties transform as 



A^A-\{aey 

i> -> i> + n e° , 



a — > a 



-+8<j>-<t>' £° 



(B21) 
(B22) 
(B23) 
(B24) 



Thus they can be unambiguously defined provided a spec- 
ified choice of time slicing, £°. Some commonly used 
gauges are |33] : the synchronous gauge in which one sets 
A = 0; the flat gauge in which -0 = 0; the longitudinal 
gauge in which er s = 0; the comoving gauge in which 
Su = 0, or equivalently Sip = in our case; the uniform 



density gauge in which Sp = 0; and the uniform Hubble 
gauge in which 9 = 3H/a, as detailed below. To com- 
pletely fix the gauge, one also needs to specify a choice 
of spatial coordinates, £. Since B and E transform as 



B ->■ B + £° - f 
E^E-Z, 
one commonly chooses to set B = 0. 



(B25) 
(B26) 



Comoving gauge: Su = Sip = 0. This gauge is de- 
fined only for <p' ^ 0. The physical degrees of freedom 
in this gauge can be represented by the gauge-invariant 
quantities 



A c = A + - (a Su)' = A + - ( a - 
a a V - 

er r = <7 S + Su = er S - 



(B27) 
(B28) 
(B29) 



The equations of motion are the simplest in this gauge, 

- \7 2 K-HV 2 <j c = -(H'-H 2 )^ A c , (B30) 
K'+nA c = 0, (B31) 
K - A c + a' c + 2U<r c = . (B32) 

Eliminating A c and cr c , one obtains a simple equation for 
1Z alone, 



K" + 2 —Tl' - c 2 V 2 TZ = 0, 

z 



(B33) 



where z = a^J '-2(H' -U 2 ) / ' c 2 s U 2 



Longitudinal gauge: er s = 0. Perturbations in this 
gauge |32j are represented by the gauge-invariant quan- 
tities 



$ = A (a a 

a 

* = V + HfJ S , 



Sy 



b'a s 



(B34) 
(B35) 
(B36) 



$ is referred to as the Newtonian potential, which is equal 
to the longitudinal curvature perturbation 4" as a result 
of Eq. (B18), $ — \& — 0. Hence it suffices to have two 



more equations of motion, 

a 2 — 

3H + H $) - V 2 * = -—8p 
*' + H$ = {H / -H 2 )Su~, 

where Su — (357), and 



(B37) 
(B38) 



8p = -- |(H'-H 2 )[ i(3t? +U8u-$) -m8v\. (B39) 
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It is easy to relate the perturbations in the longitudinal 
gauge to that in the comoving gauge. In particular, by 
Eq. (B36), we have 



(7 = Or 



and similarly, by Eqs. (B35), (B30) and (B31), 



qji _qj2 

* = K + Ucr c = - ni , V- 2 K' 
ti ci 



(B40) 



(B41) 



Synchronous gauge: A = 0. This gauge has a residual 
degree of freedom, £° = C(x)/a, where C[x) is an arbi- 
trary function of spatial coordinates. We can construct 
the following quantities that represent the perturbations 
in this gauge, 



H 
a J 


f aAdr' , 


(B42) 


a J 


f aAdr' , 


(B43) 


l_ 
a 


J aAdr' , 


(B44) 



= 5<j>- 



which are invariant up to the residual gauge freedom that 
can be absorbed into the constant of integration. The 
equations of motion for these physical degrees of freedom 
in this gauge can be expressed as 



m^ s - vV s -hv 2 <t s = --s Ps , 

= {%' -H 2 )Su s , 
4> s + <t' s + 2U(j s = , 

where 8u s = (z^z), and 

s Ps = -^(n'-n 2 )\\(5u' s + nsu s ) -msi 



(B45) 
(B46) 
(B47) 



(B48) 



By Eqs. (B42) and ( B31| , the synchronous curvature per- 
turbation tp s can be related to the comoving curvature 
perturbation 1Z through 



K + — C aA c dr' = K-- f - TZ'dr' . (B49) 
a J a J H 



Similarly, by Eqs. (B43l and (B31), the synchronous 
shear perturbation a s can be expressed in the comoving 
gauge as 



1 



Cs = Oc. — 



aAcdr' = a c + 



Uniform Hubble gauge: 5V. = -liA-tp'+\\7 2 a s = 0. 



- r^K'dr'. (B50) 
a J U 

e -HA-ij/+lV 2 <j s = 0. 
By Eq. iB6l) , the correction to the local conformal Hubble 

(B51) 



parameter at linear order is 



6n=rO-H = -UA - V' + \v 2 a s 



which transforms as 5U ->• 5% - (W -rl 2 )£° - \V 2 £° 



under the coordinate transformation (B20). Therefore, 
we can define the following gauge-invariant quantities to 
represent the perturbations in this gauge where 8% = 0, 



A H = A- 

V>H = 1p + 

a H = a s - 
84>h = S(j>- 



1 



a\{W -U 2 ) + \V 2 



5H 



H 



5H 



(w -n 2 ) + |v 2 



(w -u 2 ) + |v 2 



sn , 



sn 



(B52) 
(B53) 
(B54) 
(B55) 



The curvature perturbation Vh and shear perturbation 
<7h in this gauge can be expressed in the comoving gauge 
as 



0H 



(n'-n 2 



(w -h 2 ) + |v 2 



3 C m 



(H'-H 2 



(w -n 2 ) + \v 2 



2. vector perturbations 



(B56) 
(B57) 



There are no vector contributions to the stress-energy 
tensor from the scalar field, so the equations of motion 
are simply 



1 

Ya 
1 



2 v 2 (f; -Si) = o, 



(B58) 



(id) °(iJ)< 



which correspond to the vector modes of the (°) and (* •) 
components of the Einstein equation respectively. The 
combination 



Fa — Si 



(B60) 



is by itself gauge-invariant under the vector mode of the 
coordinate transformation 



x ->■ x + r] , 



(B61) 



where rj l i = 0. Therefore the vector shear perturbation 
aY is the one and only physical degree of freedom, which 
obeys the equations 



v V = o , 



m a) 



(B62) 
(B63) 
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